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Early SPI/INTEGRAL constraints on the morphology of the 
511 keV line emission in the 4th galactic quadrant* 
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Abstract. We provide first constraints on the morphology of the 511 keV line emission from the galactic centre 
region on basis of data taken with the spectrometer SPI on the INTEGRAL gamma-ray observatory. The data 
suggest an azimuthally symmetric galactic bulge component with FWHM of ~ 9° with a 2a uncertainty range 



covering 6° — 18°. The 511 keV line flux in the bulge component amounts to 9.9_ 21 x 10 ph cm s 



No 



evidence for a galactic disk component has been found so far; upper 2a flux limits in the range (1.4 — 3.4) x 10~ 
ph cm _2 s" 

limits on the bulge-to-disk ratio of 0.3 
Key words. Gamma rays: observations 



1 have been obtained that depend on the assumed disk morphology. These limits correspond to lower 

0.6. 



1. Introduction 

Since the first detection (Johnson & Haymes 1973) and 
the subsequent identification (Leventhal et al. 11978(1 of 
the galactic 511 keV annihilation line, the origin of galac- 
tic positrons has become a lively topic of scientific debate. 
Among the proposed source candidates figure compact ob- 
jects such as neutron stars or black holes (Lingenfelter 
& Ramaty 1983), stars, such as supernovae, novae, red 
giants and Wolf-Rayet stars, expelling radioactive nu- 
clei produced by nucleosynthesis (Ramaty, Kozlovsky & 
Lingenfelter 1979), cosmic-ray interactions with the inter- 
stellar medium (Kozlovsky, Lingenfelter & Ramaty 1987), 
pulsars (Sturrock 11971(1 . gamma-ray bursts (Lingenfelter 
& Hueter IT984J1 and stellar flares. Yet so far the source of 
the galactic positrons is still unknown. 

The question of the morphology of the galactic 511 
keV annihilation signal is intimately related to the ques- 
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tion of the origin of galactic positrons. The celestial dis- 
tributions should be tied to the source distribution, al- 
though positron diffusion within the Galaxy may to some 
extent blur this link. Although earlier measurements al- 
ready provided first indications of the morphology of the 
emission (e.g. Share et al.[1990), it is only with the advent 
of the OSSE telescope onboard the Compton Gamma-Ray 
Observatory that a first crude skymap of the 511 keV in- 
tensity distribution became available (Cheng et al. 119971 
Purcell et al. ITW71 Milne et al. EOT Milne et al. IMTTjl . 
The OSSE observations suggest at least two emission com- 
ponents, one being a spheroidal bulge and the other being 
a galactic disk component. Indications of a third com- 
ponent situated above the galactic plane have resulted 
in various speculations about the underlying source (von 
Ballmoos et al. [2003), yet the morphology and intensity 
of this component is only poorly determined (e.g. Milne 

et ai.Ennii). 

Modelling the bulge component by symmetric gaus- 
sians, the OSSE data suggest FWHM values in the range 
4° — 6° with no significant offset from the galactic cen- 
tre (Purcell et al. lTMTl Milne et al. EOT) . TGRS results 
are also consistent with no significant offset from the cen- 
tre (Ss 2°) but suggest a somewhat broader distribution 
(Harris et al. 11998(1 . The estimated bulge-to-disk ratio 
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Fig. 1. Relative exposure map of the galactic centre re- 
gion. Black corresponds to regions of maximum exposure. 

(hereafter B/D ratio) is only poorly constrained by the 
observations, and estimates vary from 0.2 — 3.3 depending 
upon whether the bulge component features a halo (which 
leads to a large B/D ratio) or not. The total flux has been 
constrained to (2.1 — 3.1) x 10~ 3 ph cm _2 s -1 (Milne et al. 
I2TTO . 

Clearly, more observations are needed to better con- 
strain the 511 keV emission morphology, and thus to shed 
light on the origin of galactic positrons. In this paper we 
report a step towards this direction by exploiting a first 
set of data recorded by the spectrometer SPI on ESA's 
INTEGRAL observatory. First results on the spectral 
shape of the 511 keV line obtained from the same set of 
data, indicating a slightly broadened line of 2.95^^ keV 
FWHM, have been reported elsewhere (Jean et al. 12003(1 . 

2. Data analysis 

The data analysed in this work were accumulated during 
the first year's Galactic Centre Deep Exposure (GCDE) 
and Galactic Plane Scan (GPS), executed as part of 
INTEGRAL'S guaranteed time observations (see Winkler 
2001J. We used data from 19 orbits from March 3rd to 
April 30th, 2003, amounting to a total effective exposure 
time of 1667 ks. The GCDE consists of rectangular point- 
ing grids covering galactic longitudes I — ±30° and lati- 
tudes b — ±10°, with reduced exposure up to b = ±20°. 
The GPS consists of pointings within the band b — ±6.4° 
along the galactic plane. For details see Winkler (2001). 
The present data are from 1199 pointings with an average 
exposure of 1400 seconds per pointing. 

As a result of data sharing agreements, the results pre- 
sented here are limited to the galactic quadrant I = 270° - 
360° but, in accordance with those agreements, data from 
pointings in the entire GCDE region I = ±30° have been 
taken into account in the analysis. The resulting sky expo- 
sure is depicted in Fig.^ where black/white corresponds 
to regions of maximum/minimum exposure, respectively. 
A quite homogeneous exposure has been achieved over 
galactic longitudes ±25°, with a small exposure dip near 
the galactic centre. The latitude dependence is approxi- 



mated by a gaussian centred on the galactic plane, with 
FWHM of - 30°. 

Data preparation and modelling of the instrumental 
background is identical to that described in Jean et al. 
(2003|). The SPI single detector event data have been 
gain corrected and binned into event spectra of 0.25 keV 
bin width for each detector and pointing, leading to a 
3-dimensional data space. The background has been mod- 
elled in this data space by two components, accounting 
for the instrumental 511 keV line and the underlying con- 
tinuum (see Eqs. 1 and 2 and Fig. 1 in Jean et al. 2003). 
The short-term variability (< 3 days) of the background 
has been predicted using the rate of saturating events in 
the Ge detectors, while the long-term variability has been 
regarded as an unknown (see section QJ. 

The expected number of counts in this data space, 
E p ,d,e, where p indicates the pointing, d the detector, and 
e the energy bin, is given by 

E p , d , e = Y,R l p^lfi + B pAe , (1) 

l,b 

where Rp b d e is the instrumental response matrix (i.e. the 
number of counts in a data space bin that arises from a 
source of unit flux at the sky position I, b; see Sturner et 
al., this volume), is the 511 keV sky intensity distribu- 
tion, and Bp d j£ is the background model. In the following 
we describe two approaches that we used to extract infor- 
mation about <&i^b from the measured SPI data. 

3. Imaging 

A qualitative impression about the spatial distribution of 
the 511 keV 7-ray line emission can be obtained by decon- 
volving the data into a celestial intensity distribution. For 
this purpose, we employed in this work the Richardson- 
Lucy algorithm (Richardson 119721 Lucy I1974fl that has 
been successfully applied to 7-ray data of earlier mis- 
sions (e.g. Knodlseder et al. 115551 Milne et al. l2HUU|l . The 
Richardson-Lucy algorithm decomposes the sky into a grid 
of equally sized pixels (here of 0.5° x 0.5°) and solves for 
the intensity "3?/^ in each of the pixels simultaneously us- 
ing an iterative maximum likelihood scheme. The pixels 
are constrained to positive intensities. As all pixel-based 
deconvolution algorithms, the Richardson-Lucy scheme 
leads to morphology artefacts in the case of reconstruct- 
ing a diffuse low significance signal since the number of 
free parameters (i.e. the number of image pixels) greatly 
exceeds the information available in the data. In order to 
reduce these artefacts, we added a smoothing step to the 
iterative scheme that effectively combines adjacent pixels 
and thus reduces the number of pixels in the reconstruc- 
tion. As smoothing kernel, a boxcar average of 6° x 6° has 
been employed. 

Figure [21 shows the resulting skymap for negative lon- 
gitudes (according to the above cited data right agree- 
ments we deliberately limit the image to I < 0°). As in- 
sets, longitude and latitude profiles of the 511 keV line 
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emission are shown that have been obtained by integrat- 
ing the intensity of latitudes b — ±5deg and longitude 
I = —5° — 0°, respectively. The emission is concentrated 
in an azimuthally symmetric region near the galactic cen- 
tre, with an extent of about ~ 9° (FWHM). The emission 
maximum is slightly offset from the galactic centre, and 
is situated at about I = — 1° and 6 = 2°. Integrating the 
intensity over the feature suggests a flux of the order of 
1CT 3 ph cm^s" 1 . 

To asses the significance of the morphology that is seen 
in the skymap, we performed Monte-Carlo simulations of 
the deconvolution process assuming that the 511 keV line 
emission morphology is described by an azimuthally sym- 
metric gaussian bulge of 10° FWHM centred at I = 0° and 
6 = 0°. In the ideal case the image deconvolution should 
reproduce the input image of an azimuthally symmetric 
emission centred at the galactic centre, yet the limited 
statistics may introduce some uncertainties in the recon- 
struction procedure. 

Figure [3] shows the resulting simulated skymap, now 
shown for the entire galactic centre region. Obviously, al- 
though the underlying model has been azimuthally sym- 
metric, the deconvolved image shows a clear emission 
asymmetry, with FWHM of 8° and 12° in the longitudi- 
nal and latitudinal directions, respectively. Also the emis- 
sion maximum is displaced from the galactic centre, and 
is found at I = 0.5° and b = 2.5°. These values are com- 
parable to those found for the centroid of the emission in 
the 511 keV sky map, and as we will show more quanti- 
tatively in the next section, the apparent displacement in 
the 511 keV sky map could indeed be of purely statistical 
nature. 



4. Model fitting 

A more quantitative approach is possible by fitting to the 
data a model of the spatial distribution that has one or 
more components. The components can be point sources, 
gaussian or other geometric forms, or arbitrary maps cor- 
responding to known distributions. The intensities of the 
components are adjusted to maximise the likelihood that 
the model gives rise to the observed distribution of counts 
in the line (a 5 keV wide energy band centred at 511 keV 
was used), binned by detector and by pointing. Along with 
the model intensities, 19 background model scaling factors 
have been adjusted for the line component by the fit, one 
for each orbit (factors G in Jean et al. I2U(J3[1 : the scaling 
factors for the continuum component (factors F in Jean 
et al. 121) (J3f) have been fixed to unity. Not fitting the back- 
ground model introduces systematic uncertainties in the 
analysis that considerable biases the morphology determi- 
nation. 

Guided by the imaging analysis (and by previous work 
performed to describe the morphology of the OSSE obser- 
vations, e.g. Purcell et al. I1997f) we used gaussian func- 
tions to describe the observed bulge emission. Assuming 
an azimuthally symmetric gaussian centred on I — 0° and 
6 = 0° results in an optimum FWHM of 9^g degrees, 
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Fig. 2. 511 keV gamma-ray line intensity map of the 
galactic centre region (only negative longitudes). Black 
corresponds to regions of maximum 511 keV line intensity. 
Longitude and latitude profiles, integrated over b = ±5° 
and I = —5° — 0°, respectively, are shown as insets. 
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Fig. 3. Simulated 511 keV gamma-ray line intensity map 
based on an azimuthally symmetric gaussian model of 10° 
FWHM centred at I = 0° and 6 = 0°. 



where the quoted uncertainties are formal statistical 2a er- 
rors. Systematic uncertainties, due to different treatment 
of the instrumental background, hardly affect the lower 
boundary, while adding an additional uncertainty of ±2° 
to the upper boundary. The flux in the bulge component 
amounts to 9.9^2 ^ x 10~ 4 ph cm _2 s _1 , where the uncer- 
tainty includes (and is dominated by) the uncertainty in 
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the gaussian width. Formally, the statistical detection sig- 
nificance of the 511 keV line amounts to 12a. 

Relaxing the condition on spherical symmetry does not 
significantly improve the fit. Yet relaxing the condition on 
the location of the bulge centroid improves the likelihood 
by 3.3 corresponding to a significance of 2.1a for the dis- 
placement. The optimum centroid positions has been de- 
termined as to be I = -1.0° ± 1.3° and b = 1.4° ± 1.3°, 
the optimum gaussian width at this location amounts to 
8^3 degrees (2a statistical errors). Systematic uncertain- 
ties in the centroid determination amount to about ±0.2°. 
The flux in the displaced bulge component amounts to 
10.312 2 x 10~ 4 ph cm _2 s _1 , including again the uncer- 
tainty in the gaussian width. 

Including a galactic disk component in addition to the 
gaussian bulge component does not significantly improve 
the fit. In none of the considered cases did the fit attribute 
a significant flux to the disk component. From the fits we 
derive upper limits on the disk flux by multiplying the 
statistical uncertainty in the disk component by the re- 
quested significance level (we quote here 2a upper limits, 
hence the formal la statistical errors have been multiplied 
by a factor of 2). Although the disk models we tested for- 
mally cover the entire galactic plane, the effective exposure 
of our data is restricted to I « ±40°, and hence we can 
only derive conclusions about this longitude range. For the 
bulge component we used an azimuthally symmetric gaus- 
sian centred on I — 0° and 6 = 0° with the width being a 
free parameter of the fit. In all considered cases, the best 
fitting FWHM amounted to ~ 9°, with a 2a lower limit 
of 6°. 

The resulting limits on the disk flux strongly depend 
on the spatial distribution that has been assumed for the 
disk component. Using models of constant positron annihi- 
lation surface density throughout the galactic disk, limited 
to a maximum galactocentric radius of 14 kpc, provide the 
largest upper flux limits, comprised between (2.8 — 3.4) x 
10~ 3 ph cm _2 s _1 for assumed exponential scale heights 
of 90 and 325 pc, respectively. A tracer of the old stellar 
population, such as the DIRBE 35 fan allsky map, gives 
a slightly smaller limit of 2.5 x 10 -3 ph cm _2 s _1 , while 
a massive star tracer, such as the DIRBE 240 /jm allsky 
map, provides a considerably smaller limit of 1.4 x 10~ 3 
ph cm^ 2 s _1 . 

5. Discussion 

The 511 keV line emission detected by SPI from the 
Galaxy is so far adequately described by a gaussian shaped 
bulge of ~ 9° FWHM. The 2a lower limit on the bulge size 
amounts to 6°, which is at the upper limit of the values 
suggested by the OSSE observations (4° - 6°). At this 
early stage of the analysis we do not emphasise this dis- 
crepancy, yet it may be taken as a first hint of morphology 
differences between the OSSE and SPI analyses. 

The data do not suggest that the bulge geometry de- 
viates from spherical symmetry yet a small offset of the 
centroid from the galactic centre direction is indicated by 



the data at the ~ 2a confidence limit. The present data 
do not yet allow to make a statement about the reality of 
the positive latitude enhancement that has been reported 
by Purcell et al. ifTMTjl . 

Assuming that the bulge is indeed located at the 
galactic centre, at a distance of 8.5 kpc, the measured 
511 keV line flux converts into an annihilation rate of 
(3.4 — 6.3) x 10 42 s _1 . The observed flux is compatible 
with previous measurements that have been obtained us- 
ing telescopes with small or moderate fields-of-view, yet 
it is on the low side when compared to OSSE measure- 
ments. OSSE, however, has detected an additional galac- 
tic disk component that is (so far?) not seen in the SPI 
data. The disk flux determined by OSSE lies in the range 
(0.7 — 2.7) x 10~ 3 ph cm _2 s _1 , and strongly depends on 
the assumed shape of the bulge component. Our 2a up- 
per limits on the disk component of (1.4 — 3.4) x 10~ 3 
ph cm _2 s _1 are still compatible with the OSSE measure- 
ments. 

The upper flux limits on the galactic disk component 
may be converted into lower limits for the bulge-to-disk 
ratio. While the constant surface density models provide 
the smallest limits of B/D > 0.3, the DIRBE 35 fim and 
240 ^m maps lead to B/D > 0.4 and B/D > 0.6, re- 
spectively. We included in these limits the 2a uncertainty 
about the bulge size which translates into a typical uncer- 
tainty of 0.2 in the bulge-to-disk ratio. Again, these limits 
are compatible with the OSSE measurement of values in 
the 0.2 — 3.3 range. 

6. Conclusions 

The analysis of SPI data that have been recorded dur- 
ing the first half of the first year's INTEGRAL GCDE 
have provided initial constraints on the morphology of the 
galactic 511 keV line emission. The data suggest that the 
emission follows an azimuthally symmetric galactic bulge 
of ~ 9° typical FWHM, yet the uncertainties on the width 
are still rather large (6° — 18°, 2a). The bulge seems cen- 
tred on the galactic centre, yet a marginal displacement 
towards positive latitudes and negative longitudes may be 
indicated. Obviously, this displacement clearly needs con- 
firmation by the analysis of a much larger set of data. 

The available SPI observations are so far rather in- 
sensitive to a galactic disk component, and only upper 
limits have been derived. More data along the galactic 
plane are needed to better constrain the disk component. 
These data are, to some extent, already taken, yet the 
data sharing agreements do not allow for their inclusion 
in the present analysis. Yet we are optimistic that once 
combined, the complete set of SPI observation will pro- 
vide unprecedented constraints on the morphology of the 
511 keV line emission, and thus give us key information 
about the origin of the positrons in the Galaxy. 
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